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Abstract 

Shielding effects in non-degenerate and degenerate plasmas are compared. A detailed derivation 
of the Wigner-Poisson system is provided for electrostatic quantum plasmas where relativistic, spin 
and collisional effects are not essential. Later a detailed derivation of a quantum hydrodynamic 
model starting from the Wigner-Poisson system is shown. The route for this derivation considers 
the eikonal decomposition of the one-body wavefunctions of the quantum statistical mixture. The 
merits and limitations of the resulting quantum hydrodynamic model are discussed. 
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I. INTRODUCTION 



Plasma physics is frequently not a part of the standard physics curriculum. Thus, theo- 
retical physicists who study quantum mechanics are not usually exposed to plasma science 
and have no interest in it. And vice versa, plasma physicists often do not have the required 
background in quantum theory. 

Plasma physics is commonly considered to be a purely classical field. However, in the 
last ten years there has been a renewed interest on plasma systems where quantum effects 
are important, see l|- 4( for extensive reviews. Quantum mechanics becomes relevant in 
plasmas when the de Broglie wavelength of the charge carriers is comparable to the inter- 
particle distance, so that there is a significant overlap of the corresponding wavef unctions. In 
such a situation, the many-body problem is described by Fermi-Dirac statistics, in contrast 
to the usual laboratory and space plasmas which obey Maxwell-Boltzmann statistics. The 
fermionic character becomes prominent for sufficiently dense plasmas. Among these, one can 
cite plasmas in compact astrophysical objects such as in white dwarfs and the atmosphere 
of neutron stars 5] or in the next generation intense laser-solid density plasma interaction 
experiments .6]. Moreover, many-body charged particle systems can not be treated by pure 
classical physics when the characteristic dimensions become comparable to the de Broglie 
wavelength. This is the case for quantum semiconductor devices, like high-electron-mobility 
transistors, resonant tunneling diodes or superlattices. The operation of these ultra-small 
devices rely on quantum tunneling of charge carriers through potential barriers. For these 
systems both Maxwell-Boltzmann or Fermi-Dirac statistics can be applied, according to the 
particle density. For instance, the_Fermi-Dirac character should be taken into account in 

I. In addition, due to recent advances in femtosecond 



the drain region of n + nn + diodes 
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pump-probe spectroscopy, quantum plasma effects are attracting attention in the physics of 
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metallic nanostructures and thin metal films 9j. Also, X-ray Thomson scattering in high 
energy density plasmas provide experimental techniques for accessing narrow bandwidth 



spectral lines 



10| , so as to detect frequency shifts due to quantum effects. Such modifications 



of the plasmon dispersion relation are inline with theoretical predictions 



precise experiments have been suggested 



111 ]. In this context, 



121 ] in order to measure low-frequency collective 



oscillations (ion-acoustic waves) in dense plasmas, as soon as keV free electron lasers will be 
available. 

The purpose of this short introduction is two- fold. First, we present an elementary 
discussion on the screening of a test charge in non-degenerate and degenerate plasmas. 
Shielding is one of the most emblematic collective effects in plasmas, so that it is worth 
to comment on the differences between the classical and quantum regimes of it. In this 
way a convenient introduction to quantum plasmas is provided. Second, we consider the 
basic kinetic and fluid models for quantum plasmas. Hence we present the Wigner-Poisson 
system, which is the quantum equivalent for the Vlasov-Poisson system in classical plasma 
physics. Afterwards we discuss the derivation of fluid equations from the Wigner-Poisson 
model. Special attention is paid to the transition from microscopic to macroscopic descrip- 
tions. While none of these subjects is new, most of the calculations here are intended to be 
shown in more detail than usually in the literature. In order to restrict the treatment to 
a basic level, only electrostatic, weakly coupled and non-relativistic quantum plasmas are 
considered. Hence, spin, magnetic field, collisional or relativistic effects are not addressed. 

This work is organized as follows. In Section II the basic physical parameters for non- 
degenerate and degenerate plasmas are introduced, by means of a discussion of the properties 
of shielding in these systems. Section III introduce the basic kinetic model for electrostatic 
quantum plasmas, namely the Wigner-Poisson system. Some of the basic properties of the 
Wigner function are analyzed, the Wigner function playing the role of a (quasi)-probability 



distribution in phase space. Section IV apply an eikonal, or Madelung decomposition of the 
one-body wavefunctions defining the quantum statistical ensemble. In this way the pressure 
functional is shown to be the sum of two terms associated to kinetic and osmotic velocity 
dispersion (to be defined later) plus a Bohm potential term associated to quantum diffraction 
effects. Section V is dedicated to final remarks. 

II. SHIELDING IN NON-DEGENERATE AND DEGENERATE PLASMAS 

Suppose a test charge q t > added to a plasma composed by an electron gas (of number 
density n(r)) and a fixed homogeneous ionic background (of number density no). Initially, 
due to the Coulomb force the electrons would have trajectories deviated toward the test 
charge. Eventually, in the equilibrium situation a stationary cloud of negative charge would 
accumulate around q t . Then, instead of q t an external observer would see an effective, smaller 
shielded charge. This is the (static) screening, or shielding effect in plasmas, which is in 
itself a manifestation of the quasi-neutrality property: due to the electric force, any excess 
charge tend to be compensated. In the stationary regime, the electrostatic field <fi = 0(r) is 
described by Poisson's equation, 



where — e < is the electron charge, Eq is vacuum's permittivity and for definiteness the test 
charge is put at the origin. For simplicity we assume the test charge to be massive enough, 
so that it can be considered at rest. 

Assuming a Maxwell-Boltzmann statistics where the quiescent electron gas is in thermo- 
dynamic equilibrium at a temperature T, one would have 




(1) 




(2) 
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where Kb is Boltzmann's constant. 

The Maxwell-Boltzmann statistics is appropriated for dilute, or non-degenerate plasmas, 
for which the degeneracy parameter x = Tp/T <C 1. Here Tp = Ep/ns is the Fermi 



temperature, defined 



13| in terms of the Fermi energy Ep, given by 



E F = ^{^n,f/\ (3) 
2m 

where h is Planck's constant over 2ir and m is the electron mass. Since electrons are fermions 
(of spin 1/2), even in the limit of zero thermodynamic temperature it is not possible to 
accommodate all of them in the ground state, due to the Pauli exclusion principle. Hence, 
excited states are filled up until the highest energy level, whose corresponding energy is 
defined as E F . 

Assuming the scalar potential to be zero before the insertion of the test charge, we can 
linearize Eq. (JT|) to obtain 

W=-=i£U -**(,). (4) 
EqKbI £o 

The radially symmetric solution to Eq. @ with appropriate boundary conditions is the 
Yukawa potential 



6 = qt P ~rl^ 
V 4- 

where 



Ad= (WY /! (6) 
V n e 2 ) 

is the (electron) Debye length. From Eq. (jHJ), an observer would measure a very small 
potential for distances larger than the Debye length, which is inline with the linearization 
procedure. Therefore \p> is a fundamental length for dilute, non-degenerate plasmas, playing 
the role of an effective range of the Coulomb interaction. The shielding effect is a collective 
effect due to a large number of electrons around qt- 
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What happens in the degenerate case? To start answering, notice that the number 
density fl2]) follows from the zeroth-order moment of the local Maxwell-Boltzmann one- 
particle equilibrium distribution function f d = / c i(r, v) given by 



/cl(r ' v) = " <2^rJ exp 



K B T 



mv 2 



(?) 



In other words, 



n(r) = J dv/ d (r,v). (8) 

On the other hand, the simplest approach for a degenerate plasma would consider an 
uniform distribution of electrons for energy smaller than the Fermi energy, and no particles 
above the Fermi level. Then Eq. ([7]) could be replaced by 

/ci(r,v) = i -| if — e0 < Ep (9) 

with /ci(r, v) = otherwise. In Eq. Qj we have the Fermi velocity Vp = (2Ep/m) 1 ^ 2 . 
Notice the energy shift due to a non-zero scalar potential, in the same manner as for the local 
Maxwell-Boltzmann equilibrium give n by Eq. (j7|). The distribution is representative of a 



zero-temperature Thomas-Fermi 14| equilibrium, displaying equal occupation probabilities 
for energies smaller than Fermi's energy, and no particles beyond. At this point temperature 
effects are disregarded. 

Equation (jSJ) can be used to find the number density for the Thomas-Fermi equilibrium 
©, yielding 

n(r) = n (l + |^) ' . (10) 
Inserting this result into Poisson's equation ([1]) and linearizing we get 

^2, 3 ™oe 2 q t xf , f . 

V <P= — — -0- — 5(r). (11) 

which is the same as Eq. 01]) with the replacement k^T — > Ep, except for a numerical factor. 



Therefore, a shielding distance Xp, or Thomas- Fermi length 



3], can be set for degenerate 



plasmas, 

1/2 

A 



The Thomas-Fermi length is non-zero even for zero thermodynamic temperature, unlike 
the Debye length Ad- This happens because of the exclusion principle which prevents the 
accumulation of electrons in the same place as the test charge. 

Notice that in a sense every particle in a plasma, be it degenerate or not, can be in- 
terpreted as a test charge with the corresponding screening cloud. Hence instead of the 
long-range Coulomb field we have an effective Yukawa interaction field. This point of view 
is adopted, for instance, in the treatment of the ultrafast phase-space dynamics of ultracold, 
neutral plasmas 15 1. 

In our simplified picture / c i(r, v) was regarded as a purely classical probability distribution 
function. A more detailed treatment taking into account quantum diffraction effects shows 
that at large distance the inter-particle oscillation behave as 

0~^cos(2r/A F ), (13) 



a phenomenon known as Friedel oscillations ly, |17|]. As discussed in the next Section, the 



Wigner function provides a convenient tool to incorporate quantum effects in plasmas, in 
strict analogy with the classical probability distribution approach. 

The analysis of shielding allows the introduction of characteristic length scales A^> and 
\p in non-degenerate and degenerate plasmas, respectively. In addition a pertinent length 
scale measuring the spatial extension of the wavefunction of each electron is the de Broglie 
wavelength A# = h/(mvT), where vt = (2kbT '/m) 1 / 2 is the thermal velocity. Given the de 
Broglie wavelength one can express the degeneracy parameter as 

X = ^ = \^ 2 n,\\)^. (14) 



Hence Fermi-Dirac statistics is necessary when Xb is of the same order of the inter-particle 
distance. 

Other fundamental scales in classical and quantum plasmas are as follows. 
• Time scale for both classical and quantum plasmas: w" 1 , where 

2\ 1/2 



"* = \ik) • (15) 

is the plasma frequency. Given some electron charge depletion, an electric force appears 
in order to restore complete charge neutrality. The resulting linear oscillations have a 
frequency lo p . 

Typical interaction energy U- in t for both classical and quantum plasmas: 

-1/3 



U mt = 6 -^, (16) 



since the mean inter-particle distance scales as tiq 

• Typical kinetic energies: Kq = k>bT for non-degenerate and Kq = kbTf for degene- 
rate plasmas. This happens for fermions because of the filling up of excited states due 
to the Pauli exclusion principle, even at zero temperature. 

From the typical energy scales we can form classical Tq and quantum Tq energy coupling 
parameters, 

77 2 V 3 1/3 

T Q = ^ = — — — — —pr = 5.0 x 10 10 no 1/3 • (18) 

" Kq (3vr2)2/3 £o ^ n J/3 



Numerical values are for S. I. units. Weakly coupled plasmas have small energy coupling 
parameters. From Eqs. ffTTHTS]) the conclusion is that while classical weakly coupled plasmas 
tend to be dilute and cold, quantum weakly coupled plasmas tend to be dense. For example, 
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with tiq > 10 35 m~ 3 (white dwarf) one has Tq < 0.1, allowing the use of collisionless models 
in a first approximation. This is a consequence of the Pauli exclusion principle, which 
contribute to forbids e-e collisions in very dense systems. In other words, the more dense a 
degenerate plasma is, the more it resembles an ideal gas, except for the mean field, collective 
interaction. 

A more detailed account on the characteristic scales in classical and quantum Coulomb 
systems can be found 



xount on tne 
in Q,Q,Q. 



III. OBTAINING THE WIGNER-POISSON SYSTEM 

The Vlasov-Poisson system 

^ + w ^ + JL^^£i = o (19) 
dt dx mdx dv 

|^ = y {J -no) (20) 

is the basic tool for kinetic theory of classical plasmas, where for simplicity we take an 
one-dimensional electron plasma in a fixed neutralizing ionic background no- In the Vlasov 
equation f|T9|) . f c i = f c i(x,v,t) is the reduced one-particle probability distribution function. 
For a plasma with N electrons, (1/N)f c i(x, v , t) dv dx gives the probability of finding one 
electron with position between x and x + dx and velocity between v and v + dv, at the time 
t. Hence, the normalization 

J dvdxf cl (x,v,t) = N (21) 

is assumed. 

From the knowledge of f c i, obtained solving either analytically or numerically the Vlasov- 
Poisson system subject to appropriate boundary conditions, the whole machinery of classical 
statistical mechanics can be used to compute the expectation values of macroscopic quanti- 
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ties. For example, the average kinetic energy of one electron follows from 

mv 2 1 f , , „ , . mv 2 , . 

< — >= -ft / dvdxf d (x,v,t)—. (22) 

It is convenient to adopt a similar methodology in a quantum kinetic theory for plasmas, as 
far as possible. 

Quantum mechanics can be formulated in phase-space using the Wigner function / = 
f(x, v, t), which provides a quantum equivalent of f c i(x, v, t). Let us study some of the basic 
properties of the Wigner function. For an one-particle pure state quantum system with 
wavefunction ip(x,t), the Wigner function is defined 



|ia| by 

'^/M^W^'M*-?')- (23) 

with all symbols as before. From f(x,v,t) we can compute the probability density 

dvf(x,v,t) = ty(x,t)\ 2 (24) 
and the probability current 

dvf(x , v , t)v = ilL^l- r m. (25) 
2m \ ox ox J 

The above wavefunction is normalized to unity. 

In general, even an one-particle quantum system can not be represented by a wave- 
function alone. Mixed quantum states are described by a quantum statistical ensemble 
{i/j a (x, t) ,p a }, a = 1,2, ...M, with each wavefunction ip a (x,t) having an occupation pro- 
bability p a such that p a > ,J2a=iPa = 1- The Wigner function is then defined by the 
superposition 

m / / tmvs\ . ( s \ , / s \ . . 

f = ^2^P° J dsex v(^T M X+ 2'VM X -2'V ' (26) 

Q=l J V 7 
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Correspondingly, we have the probability and current densities 

/M 
dvf(x,v,t) = J2Pa\Mx,t)\\ (27) 
a=l 

The density matrix p(x, y, t) could also be elected as the central object in a quantum 
kinetic theory for plasmas. However, in this manner the similarity to the Vlasov-Poisson 
model would be lost. Besides, we have a complete correspondence since 

M 

p(x,y,t) = Paip a (x,t)ip* a (y,t) 

a=l 

f , ( imv(x -y)\ .fx + y \ 
= dvexpi ) H^~' v,t ) ' ^ 



with inverse given by 



m f fimvs\ ( s s \ 



(30) 



In other words, given the Wigner function the density matrix can be found and vice-versa. 
The same apply to many-body Wigner functions and density matrices. 

Going one step further, consider now a N— particle statistical mixture described by 
the set {ipa (x\,X2, ...,XN,t) , p a }, where the normalized iV— particle ensemble wavefunc- 
tions ip^(xi, x 2 , xn, t) are distributed with probabilities p a , a = 1,2, ...,M satisfying 
Pa > 0, J2a=i Pa = 1 as before. Here Xi represents the position of the ith-particle, 
% = 1,2, ...TV. All particles have the same mass m. The N— particle Wigner function is 
defined by 

m L,=l ^ 



/ j \ n ^ r 

f N (x 1 ,v 1 ,...,x N ,v N ,t) = N -j £^Pa ds!...ds N exp 



h 

x €*^i + S i,--^N+ S f,t)^(x 1 -^,...,x N - S f,t). (31) 
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The factor N in Eq. (I3TT) is inserted so that 



/ 



dxidv i...dxNdvN f vi, xjy, ujv, t) = N. 



(32) 



In this manner, the integral of / over all the velocities gives a number density in configura- 
tion space. The fermionic character of the system is not necessarily included, but this could 
be done assuming the N— body ensemble wavefunctions to be antisymmetric in Eq. (l3Tj) . 

Wigner functions provide a convenient mathematical tool to calculate average quantities, 
in the same way as for the classical probability distribution function. However, they are 
not necessarily positive definite. Hence it is customary to refer to them as quasi-probability 
distributions. 

To proceed to the derivation of the quantum analog of the Vlasov-Poisson system (fTDH^Ul) , 
we can introduce the reduced one-particle Wigner function f(xi,v\,t), 



and the reduced two-particle Wigner function f^ 2 \xi, v i, x%, t>2, t) with a convenient norma- 
lization factor N, 



Similar reduced N— particle Wigner functions with N > 3 can be likewise defined. If the 
Wigner function were a true probability distribution, (1/iV) f(xi, v±, t)dx\dv\ would give 
the probability of finding the particle 1 in an area dx\dv\ centered at (xi,t>i), irrespective 
of the "position" and "velocity" of the remaining ith-particles, % = 2,...,N. An analogous 
"classical" probability interpretation could be assigned to the remaining reduced Wigner 
functions, except for the non-positive-definiteness property. 
In passing, we have 




(33) 



f {2) (x 1 ,v 1 ,X2,v 2 ,t) = N / dx 3 dv 3 ...dx N dv N f N (x 1 ,v 1 ,...,x N ,v N ,t) . (34) 




(35) 
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It is reasonable to pay special attention to the reduced Wigner functions, since f N contain 
far more information than what is commonly needed. In this regard the one-particle Wigner 
function plays a distinguished role. Indeed, macroscopic objects like number and current 
densities can be derived from / after integration over just one velocity variable, exactly as 
in Eqs. fl27H28|) . originally written for an one-particle system. 

re evolution equation satisfied by the one-body Wigner function, the philo- 



To obtain t 
sophy of Ref. 



20[ can be pursued. Consider then the Schrodinger equation satisfied by the 



iV— body ensemble wavef unctions, 



ih- 



Ot 



2m ^— ' dx 

i=i 



d 2^N 

Y + V(xi, ...,x N ) V> 



(36) 



for an interaction energy V(xi, ...,Xn). 

We are concerned with the case where the system components interact through some 
two-body potential W, 

V(x 1 ,...,x N ) = Y,W(\x l -x ] \). (37) 

i<j 

This situation is evidently interesting because of the Coulomb force. 



Some algebra shows 



d£ + df_ 
dt dxi 



20| that 



i m 



ds\ dv[ dx<i dv' 2 exp 



x (W{\x x -x 2 + ||) - W(\ Xl -x 2 - ||)) f^(x 1 ,v[,x 2 ,v' 2 ,t) 



(38) 



in terms of the reduced two-particle Wigner function f( 2 \ In the derivation, iV 3> 1 
was taken into account. Actually a more detailed argument involving the higher-order 
Wigner functions yield a quantum BBGKY (Bogoliubov-Born-Green-Kirkwood-Yvon) hie- 



rarchy 



21 



24|, where the dynamics of the (N — l)-body reduced Wigner function is shown 
to depend on the iV-body reduced Wigner function. Hence in both the classical infinite 
BBGKY set of equations and its quantum analogue we are faced with a closure problem. 
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Ignoring correlations is the simplest way to close the system, considering that the distribu- 
tion of particles at (xi, Vi) is not affected by particles at a distinct phase space point (xj, Vj). 
In this mean field (or Hartree) approximation the N— body Wigner function factorizes, 

f {2) (x 1 ,v 1 ,X2,v 2 ,t) = f(x 1 ,vi,t) f(x 2 ,v 2 ,t) . (39) 



Equation ([35]) then becomes 

d f d f f 

Jt +Vl dx^ = I dv'iKWscWx-VnX^fix^t), (40) 
with the mean field self-consistent potential 

W sc (x,t) = J dvdx' f(x',v,t)W(\x- x'\) . (41) 

The functional i^[Wg C | v[ — Vi, Xx, t] is defined by 

T ^ rTT , , / n im f , ( im(v\ — V\)s\\ 
KlWsclv^-v^xut) = -^-fp / dsi exp I Ix 

x (W SC ( X1 + j,t)- W sc (xt - ^, t)) . (42) 

Frequently an external, possibly time- dependent potential V ext (xi, ...,XN,t) should be 
included. For instance, such a circumstance arises in solid state devices, when considering 
the electronic motion in a fixed ionic lattice or under a confining field like in quantum wires 



or quantum wells 



a la 



251 ]. Or even the field due to an homogeneous ionic background can 
be thought as an external superimposed field. Hence, consider an external potential of the 
form 

N 

V ext {x u x N , t) = ^ W ext {xi, t) (43) 

i=l 

for some one-particle potential W ex t(xi,t). Implicitly in Eq. (T43|) . the functional form of 
W ex t is the same irrespective of X{, implying that the external field has the same influence 
on all particles. For completeness we indicate the changes for a potential 

N 

V{x 1 ,...,x N ) = W{\xi - + Y, W ext {xi,t) . (44) 

i<j i=l 
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Following the same steps as before, the one-body reduced Wigner function f(xi,Vi,t) can 
then be shown to satisfy 

Jt +Vl ld~ 1 = I dv i K l W ^ + W ext \v' 1 -v 1 ,x 1 ,t}f{x 1 ,v[,t), (45) 

for 

T^r-rrr TT , , / n im f , ( im(v\ — t>i)si\ 
K[W SC + W ext \v 1 -v 1 ,x 1 ,t] = -^fc2 / dsi ex P ( ^ ) x 

x (V sc (xi + y , t) + W ext (xi + j,t)- W sc (xi - t) - W ext ( Xl - |-, t)) 

(46) 

and for the averaged self-consistent potential W sc the same as in Eq. ( ]4"TT) . 

The necessary changes in three-dimensional charged particle motion are as follows. As- 
sume the Coulomb interaction 

W(\r-v'\)= f (47) 
47T5 |r — r'| 

and define the total electrostatic potential 4>(r,t) so that 

<f>(T,t) = <f> BC (r,t) + <l> ex t(T,t), (48) 
in terms of the self-consistent W sc and some external W ext potentials, where 

W sc {r, t) = -e0 sc (r, t) , W ext (r, t) = -e<p ext (r, t) . (49) 
From the three-dimensional version of Eq. fj4~Tj) it follows that 



_± f dvdv'f(r',v,t)\7 2 ( p 
£ J V4vr|r-r'| 

— / dv dr' f(r',v,t) 8(r — r') 
£o J 

J dv/(r,v,f). (50) 



e 
s 



15 



Moreover, 



V 2 ^ = ~-V 2 W ext = -™ (51) 

e e 

ixed homogeneous ionic background of 



25] are needed in the case of a non- 



if the external potential is due to an immobile 
density n and ion charge e. Appropriate changes [7|, |8|, 
homogeneous background e.g. as in the case of doped semiconductors, or in the presence of 
a dispersive medium with a permittivity constant e 7^ Eq. 
From Eqs. floUHBT]) it is immediate to derive 

V ^ = io (/ dY ^ r ^^- n °) ' ( 52 ) 
the Poisson equation in this case. 

For notational simplicity, it is indicated to restrict again to the one-dimensional case. In 
terms of the electrostatic potential 4>, Eq. (J4"5j) is rephrased as 

dt dx 



df f 

+ V ~dx = dv K<f,[<j)\v -v,x,t]f(x : v',t) , (53) 



where K^> \ v' — v, x, t] is the functional 



_ r , 1 / n iem f ds ( im(v' — v)s\ 
K^[(p \v — v , x, t\ = — — / — - exp I x 



h J 2nh \ h 
x fax+irf-ftx-ltj) . (54) 

Equation ( 1531) can be termed the quantum Vlasov equation (in the electrostatic case), since it 
is the quantum analog of the Vlasov equation satisfied by the reduced one-particle probability 
distribution function. Finally, the quantum Vlasov equation should be coupled to Poisson's 
equation, 

dv f(x, v, t) — no ) • (55) 



8x 2 Eq 

Equations ( 15 3p and fl5lj|) constitute the Wigner- Poisson system, which is the fundamental 
kinetic model for electrostatic quantum plasmas. It determines in a self-consistent way both 
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the Wigner function, associated to how the particles distribute in phase space, and the scalar 
potential, which in turn describe the forces acting on the particles. 

The Wigner- Poisson system needs to be supplemented with suitable boundary and initial 
conditions. For plasmas, decaying or periodic boundary conditions are frequently employed. 
For nano-devices, the choice of boundary conditions is subtler due to the finite size of the 
system and the nonlocal character of the Wigner function. Indeed, to compute the integral 
defining the Wigner function we need to specify f(x,v,0) in the whole space even when 



dealing with finite size systems 



7. 8 
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Let us review the necessary steps for the derivation of the Wigner-Poisson system. Above 
all, it is a mean field model with the N— body ensemble Wigner function supposed to be 
factorisable, in order to achieve the simplest closure of the quantum BBGKY hierarchy. 
Hence the N— body Schrodinger equation (or the Liouville-von Neumann equation for the 
N— body ensemble density matrix) is replaced by a system with fewer degrees of freedom. 
Indeed, for N electrons in three-dimensional space, we have 3N + 1 coordinates to define 
the wavefunction, 6iV + 1 coordinates for the density matrix, and only 6 + 1 = 7 indepen- 
dent variables for the reduced one-body Wigner function, here taking into account time. 
Therefore, the mean field theory is much less numerically demanding, since it deserve the 
discretization of a space with fewer dimensions. The price for the reduction is the neglect 
of collisions, besides spin and relativistic effects at least in the present formulation. Finally, 
no magnetic fields were included. 

Since it is the analog to the Vlasov- Poisson, the Wigner-Poisson system becomes the natu- 
ral tool in quantum kinetic theory for electrostatic plasmas. This is because the methods ap- 
plied to the Vlasov- Poisson system can with some optimism be translated to Wigner-Poisson 
quantum plasmas. Nevertheless, other quantum kinetic treatments for charged particle sys- 



tems are obviously important. For example, the density functional 

17 



251 ] and Green's function 



26l . l27j ] approaches are popular tools for the modeling of quantum transport among the 
condensed matter community. Moreover, the simplifications of the Wigner-Poisson model 
can be sometimes overcomed by alternative formulations. For instance, Green's function 
techniques can be used to describe collisions associated to short range particle-particle in- 



teractions 
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27] , in terms of a Boltzmann type collision operator. 
It is interesting to look to the semiclassical limit of the quantum Vlasov equation (15H|) . 
By means of the change of variable s = hr/m and Taylor expanding, one get 

df df e d(f)df = eh 2 <9 3 </> <9 3 / Q 4 
dt dx mdxdv 24 m 3 dx 3 dv 3 

Implicitly, the semiclassical approximation assumes the smallness of a non-dimensional quan- 
tum parameter H = fij (mvoLo), where Vo and Lq are resp. characteristic velocity and length 
scales. 

Equation ( )56|) is a semiclassical Vlasov equation, with / playing the role of one-particle 
distribution function. We see that unlike for classical plasmas, in general neither / nor phase 
space volume are preserved by the quantum Vlasov equation, since 

«**±*L + 0W + S (57) 
dt 24m 3 dx 3 dv 3 

along the (classical) characteristic equations 

dx dv e d(f> , , 

dt ' dt m dx 

Due to this property the positive definiteness of the Wigner function is not preserved by 
Eq. ( )53l) . except for linear electric fields and a vanishing quantum correction. Also notice 
that Eq. (157)1 is a Boltzmann's like equation, although this is not exactly true since one has 
time-reversal invariance under t — > —t,x — > x,v — > —v, f(x,v,t) — > f(x, —v, —t) (assuming 
4>(x,t) = (p(x,—t)). There is no irreversibility nor memory loss at all in the quantum 
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Vlasov equation, be it in the semiclassical or the fully quantum versions. This is not exactly 
surprising since the Schrodinger equation is time- reversal invariant. 

Even when the quantum Vlasov and Vlasov equations are the same, which happens for 
linear electric fields, f(x,v,t) can not be considered as an ordinary probability distribution 
function. Not all functions on phase space can be taken as Wigner functions, since a genuine 
Wigner function necessarily correspond to a positive definite density matrix. Hence, we have 



at least 



281 ] the following necessary conditions, 



/ 



dxdvf = N , (59) 

dvf > 0, (60) 

dxf > 0, (61) 

< (62) 

Equation (|59|) is simply a normalization condition, while Eqs. (!60|) and ( l6Tj) assure the 
spatial and velocity marginal probability densities to be everywhere non-negative. Equation 
(l62p eliminate too spiky Wigner functions, which would be against the uncertainty principle. 
For instance, for the Gaussian profile 

' = ^ exp HJ) exp (~S) (63) 

with constant standard deviations a x>v from Eq. (l6"2~j) it follows that 

ft 

CTa:^ > 7T~ ) ( 64 ) 

in accordance with the uncertainty principle. 

Since no collisional effects are included in the Wigner- Poisson system, strongly coupled 
quantum plasmas deserves special treatment. Hence in principle the Wigner- Poisson model 
assumes a small energy coupling parameter, see Eqs. fjTTHT8|) . However, for dense plasmas 
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(as the electron gas in metals) sometimes quantum collisionless models are still applicable 
thanks to the Pauli blocking phenomenon preventing e-e collisions {4], even if the energy 
coupling parameter is large. 

Not only very dense charged particle systems deserve quantum kinetic equations. For 
instance, due to the ongoing miniaturization, even scarcely populated electronic systems such 
as resonant tunneling diodes 7j should be described in terms of quantum models. Indeed, 
the behavior of these ultra-small electronic devices relies on quantum diffraction effects as 
tunneling, making purely classical methods inappropriate. The non-local integro-differential 
potential term in Eq. ( 15 3 p in the Wigner-Poisson system has been shown to be capable of 



the modeling of negative differential resistance, associated to tunneling 29]. Moreover, the 
collisionless approximation becomes more reasonable in view of the nanometric scale of the 
devices, simply because the mean free-path exceeds the system size. In the same manner, 
the usually extreme high operating frequencies makes the collisionless approximation more 
accurate, because ojt <C 1 for an operating frequency 00 and a average time r between 
collisions. For example, in resonant tunneling diodes one can find [3] potential barriers 
of the order 0.3 eV ~ hw, implying an operating frequency u ~ 10 15 s _1 . Therefore the 
Wigner-Poisson system is well suited for ballistic, collisionless processes in nanometric solid 
state devices, even at relatively low densities no ~ 10 24 m~ 3 . Correspondingly one finds a 
Fermi temperature Tp ~ 40K much smaller than a typical room temperature T ~ 300i^, 
justifying the non-degeneracy assumption and Maxwell-Boltzmann's statistics. 

IV. DERIVATION OF A FLUID MODEL FOR QUANTUM PLASMAS 

The Wigner-Poisson method presents some drawbacks: (a) it is a nonlocal, integro- 
differential system; (b) its numerical treatment requires the discretization of the whole phase 
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space. Moreover, as is often the case with kinetic models, the Wigner-Poisson system gives 
more information than one is really interested in. 

For these reasons, it would be useful to obtain an accurate reduced model which, though 
not providing the same detailed information as the kinetic Wigner-Poisson, could still be 
able to reproduce the main characteristics of quantum plasmas. 

To obtain a set of macroscopic equations for quantum plasmas, first we derive a system 
of reduced 'fluid' equations by taking moments of the Wigner-Poisson system. Using a 
Madelung (or eikonal) decomposition, it can be shown that the pressure term appearing 
in the fluid equations can be separated into a classical and a quantum part. A working 
hypothesis is then applied to the classical term, so as to close the fluid system. Shortly the 
meaning of classical and quantum contributions to the pressure will be explained. 

This approach to a quantum hydrodynamic model for plasmas appeared in Ref . 30( . The 



same quantum fluid model has been applied to several distinct problems involving charged 

n 

particle systems, for instance, the nonlinear electron dynamics in thin metal films [9(, the 



3l|, parametric amplification characte- 



excitation of electrostatic wake fields in nanowires 

ristics in piezoelectric semiconductors [32;], breather waves in semiconductor quantum wells 
33] ]. multidimensional dissipation-based Schrodinger models from quantum Fokker-Planck 
dynamics 



34|, the description of quantum diodes in 



degenerate plasmas 



35| and quantum 



to name but a few. The extension 



33| 



ion-acoustic waves in single- walled carbon nanotubes 
of the model to incorporate magnetic fields was done in 

We take moments of Eq. (1531) by integrating over velocity space. In other words, intro- 
ducing the standard definitions of density, mean velocity and pressure 

n(x, t) = I f dv , u(x, t) — — I fv dv , P(x, t) — m ( ( fv 2 dv — nu 2 j , (65) 
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we get 

dn d(nu) , . 

du du e 86 1 dP 

^T + u lT = -7r 7T- ( 67 ) 

at ox max mn o x 

A more detailed theory would include the energy transport equation obtained after taking 
the second-order moment of the Wigner function and the associated time-derivative. 

Equations (I66T4671) do not differ from the ordinary evolution equations for a classical 
fluid. This may seem strange, but in the following it will appear that the quantum nature of 
the system is in fact hidden in the pressure term. Contributions where h explicitly appear 
can be found only in the higher-order moments. Actually, taking into account the energy 
transport equation is not sufficient, since in the electrostatic case h appear explicitly only 
in the equation of motion for the third-order moment 3f|. 

To proceed, first notice the equivalence between the Wigner-Poisson and a system of 
countably many Schrodinger equations coupled to the Poisson equation, as has been ma- 
thematically demonstrated [39|. More exactly, the reduced one-body Wigner function can 
always be written as 

f( x,v,t ^ = 2~n 5^ Pa I ds 6Xp \ ~~T~ ) + ~ |,*) > ( 68 ) 

with ensemble probabilities p a > so that J2a=iPa = 1, for each one-particle ensemble 
wavefunctions ijj a (x,t) normalized to unity and satisfying 

dip d^ip 

<^ = -_^- c ^ a> a = l,...,M, (69) 

which is the Schrodinger equation for a particle under the action of the mean field electro- 
static potential 4>(x,t). In addition, the Poisson equation flo^]) is rewritten as 

Q2 - / 



dx 2 e 



Nj^PalM^t)] 2 -n . (70) 



a=l 
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Equations (169H7D1) constitute the so-called Schrodinger-Poisson system, which has to be 
supplemented with suitable initial and boundary conditions. It provides a way of replacing 
the original iV— body problem by a collection of one-body Schrodinger equations, coupled by 
Poisson's equation. From a methodological point of view the Schrodinger-Poisson modeling 
corresponds to put the emphasis again on the wavefunction and not on the (phase space) 
Wigner function. Collective effects are mediated by the self-consistent potential <fi. 

A rigorous proof |39| of the equivalence of Eqs. ( 169H70|) and the Wigner- Poisson system is 
beyond the present text. However, at least we can obtain some insight on the interpretation 
of the ensemble wavefunctions. From Eqs. (jMHSS])? 

f(x\,Vi,t) = J dX2<iv2---dXNdVN f N (x\, V i, Xn, V TV, t) 

M , / • \ 

.. / m / , , , ( imv\8\\ 
= \2nh) 2-^ Pa dsidx 2 ...dx N exjp I - — 1 x 

a=l J ^ ' 

x V>« * \Xl + ~£> x 2> -;BN,tj ( x l ~ T> X2 > ■■■> X N, t ) ■ ( 71 ) 

To be inline with the mean field approximation it is natural to factorize as 

Ipa i X li X 2, XN, t) = 1p a (x 1 ,t) X ... X Tp a (x N , t) , (72) 

for the A^— body wavefunction, fully neglecting correlations. Quantum statistics effects are 
not taken into account in the Ansatz (1721 . which does not respect the Pauli principle. With 
this proviso, inserting Eq. (172"]) into Eq. (ITTj) the result is precisely Eq. ([55]) . with the same 
statistical weights p a . Hence we can view the one-body ensemble wavefunctions ip a (x,t) 
as the result of splitting the N— body ensemble wavefunction into the product of identical 
factors. Actually it not so surprising that a correlationless model could at the end be written 
in terms of a collection of one-body Schrodinger equations. 

Thanks to the Schrodinger-Poisson form, we are able to decompose the pressure term in 
a classical and a quantum part, as follows. Consider the Wigner distribution in Eq. (I68l) . 
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In terms of the ensemble wavef unctions, from Eq. ( 165]) one obtains 



M 



n 



2 



a=l 



~ - /.Pa\Ya a _ y a a _ ) 



2m 



a=l 



dx 



and, after some work, 
P 



Nh 2 

4m 



a=l 



dip c 



dx 



N 2 h 2 

4mn 



Now we Madelung 



M 



a=l 



d 2 r a 

dx 2 



dx 



40| decompose the wavefunctions with 



ip a (x,t) = A a (x,t) exp (iS a (x,t)/h) 



where the amplitudes A a and phases S a are real functions. We get 

M 



n 



l 2 



<*=1 



nw = — > P Q ^4 
m z — ' 

a=l 



So; 



and also 



P 



N 2 
2mn 



a,/3=l 

Nh 2 ^ 
+ —l^P* 



dS n dS, 



2m 



a=l 



d A, 



\ dx dx 

2 



A r 



d 2 A r . 



(73) 
(74) 



(75) 



(76) 



(77) 
(78) 



(79) 



dx J dx 2 

In the pressure there is now the explicit presence of h. However, notice that since the 
ensemble wavefunctions satisfy the one-body Schrddinger equation ( 1691) both the amplitudes 
and phases implicitly depend on Planck's constant. 



It is useful to define the kinetic u a and osmotic u° a velocities associated to the wavefunction 



lJ>o 



1 dS, 



u. 



m dx 



h dA a /dx 

m A n 



(80) 
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Then it can be verified that the pressure in Eq. (179|) is given by 

P = pk + po + pQ ; ( gl ) 

where the kinetic pressure P k is 

M 

pk= rrm p a pp(u a - Up) 2 , (82) 

o,/3=l 



the osmotic pressure P° is 



M 

P ° = ^ E P«P/»(««"^) a , (83) 



and the quantum pressure is 

In Eqs. (I82H831) . there is a modified set of ensemble probabilities p a = p a (x,t), 

ft. = ^. (85) 
n 

The new statistical weights satisfy p a > 0, 5^ a =i Pa = 1 as they should. 
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42|, 



For a particular ^ a the osmotic velocity points to the regions of higher density 
as is more evident in the three-dimensional version, 

u° = (h/m)V\nA a . (86) 

Both pressures P k and P° can be viewed as a measure of the dispersion of velocities, 
kinetic and osmotic. Consider the following redefined average < f a > of an ensemble 
function = f a (x,t): 

M 

< fa>=^P*f a - (87) 
a=l 

Equations f[8"2"H8"51) are equivalent to the standard deviations 

P k = mn((u 2 a ) - (u a ) 2 ) , (88) 
P° = mn(([<] 2 )-«) 2 ). (89) 
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For a pure state so that p a = 5 a/ 3 for some (3, both P k and P° vanishes and only P® survives. 

Since the kinetic and osmotic pressures are a measure of the kinetic and osmotic velocities 
dispersion, it is reasonable, although not rigorous, to assume an equation of state so that 



For defmiteness, we call P c the "classical" part of the pressure, in the sense that it represents 
a measure of the velocities dispersion. However, it explicitly contains Planck's constant since 
it depends on the osmotic velocities, which always have a purely quantum nature. 

The replacement of the sum of the kinetic and osmotic pressures by a function of the den- 
sity only requires some comments. Equation (IHTj) is exact but offer no advancement over the 
Wigner-Poisson formulation, because ultimately it requires the knowledge of the ensemble 
wavefunctions and then the solution of a countable set of self-consistent Schrodinger equa- 
tions. In classical kinetic theory (therefore, without the osmotic and quantum pressures), it 
is customary to assume a closure assuming that the standard deviation of the velocities is 
a function of density only. We have just gone one step further, including also the standard 
deviation of the osmotic velocities. 

In addition, in the classical limit we expect equations reproducing the classical fluid 
equations. This is certainly true if P k + P° = P c '(n) for some appropriate function of 
density only. Finally, the standard Euler equations are reproduced thanks to the residual 
classical limit in P k . Or, for h = 0, we also have P° = and we expect P k to be some 
function of the density only. 

Unlike P in Eq. flo^]) which uses the Wigner function, in Eqs. (jHEHEH]) the statistical 
weights are provided by the p a in Eq. fl85l) . For a pure state one has P c = 0, which is inline 




(90) 



depending only on density. In this way, we obtain 




(91) 
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with the understanding that a pure state corresponds to a cold plasma with no dispersion 
of velocities. The contribution P®, on the other hand, have no classical counterpart at all. 

It is well-known that the closure problem is a delicate one. The derivation of macroscopic 
models from microscopic models always deserve some degree of approximation and a more or 
less phenomenological or ingenuous point of view. The present approach is capable of taking 
into account the quantum statistics of the charge carriers, represented by an appropriated 
equation of state. For example, assuming interaction with a heat bath an isothermal equation 
of state is indicated. For fast phenomena where thermal relaxation has no time to occur, 
an adiabatic equation of state can be used. In addition, Maxwell-Boltzmann or Fermi- 
Dirac distributions can be applied according resp. to the dilute or dense class of plasma. 
Moreover, the model takes into account quantum diffraction effects, in particular tunneling 
and wave packet dispersion, present in the quantum part P Q of the pressure. It is also able 
to reproduce the linear dispersion relation from kinetic theory if the equation of state is 
adequate, except for purely kinetic phenomena. Finally, the quantum fluid model reduces 
to the standard Euler equations in the formal classical limit and are sufficiently simple to 
be amenable to efficient numerical simulation. 

The proposed simplification becomes more justified for an important class of statistical 
ensembles, where the wavefunctions have all equal (but not necessarily constant) amplitude, 



Then the osmotic pressure identically vanishes and one has p a — p a , so that Eq. flHHj) 
becomes the usual standard deviation of the kinetic velocities. Hence P k can be interpreted 
in full analogy with the standard thermodynamic pressure. The velocities dispersion arises 
just from the randomness of the phases of the wavefunctions. The approximation can be 
viewed as a first step beyond the standard homogeneous equilibrium of a fermion gas, for 




(92) 
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which each state can be represented by a plane wave 

ip a (x, t) = A Q exp (imiiax/h) , (93) 

with the amplitude A and the velocities u a spatially constant. In the generalization (J92]), 
both the amplitude and the velocity can be spatially modulated, although the amplitude 
is the same for all states. For systems not so far from equilibrium, this appears to be 
reasonable. 

With the hypothesis f[9"Uj) the force equation (16T|) can be written as 

du du 1 dP c (n) e d (f> 1 d P Q 
at ox mn ox max mn d x 



Using the identity 

1 dP Q h 2 d fd 2 {^i)/dx 2 



(95) 



mn d x 2m 2 dx \ \fn 
we can rewrite the basic quantum hydrodynamic model for plasmas in terms of the continuity 
equation fl66l) and the force equation 

du +u dv L 1 dP c {n) e_ ft 2 d f d 2 {^n)/dx 2 \ 
dt dx mn dx mdx 2m 2 dx \ ^fn ) 

In the limit h — > this is formally equal to Euler's equation for an electron fluid in the 
presence of an electric field —d<ft/dx. Finally, we have the Poisson equation 

° 2 " 6 in -n ), (97) 



dx 2 e 

The only difference to classical plasmas is the ~ h 2 term in Eq. ( )96l) . the so-called Bohm 
potential term. While mathematically the Bohm potential is equivalent to a pressure to be 
inserted in the momentum transport equation, physically it corresponds to typical quantum 
phenomena like tunneling and wave packet spreading. Therefore it is not a pressure in the 
thermodynamic sense. Besides, it survives even for a one-particle pure state system. 
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Many quantum hydrodynamical models are popular in the context of semiconductor 



physics. For instance Gardner 
introduced by Wigner 



43| considered a quantum corrected displaced Maxwellian as 
19j . More exactly, it is possible to find the leading quantum correction 
fi(x,v,t) for a momentum-shifted local Maxwell-Boltzmann equilibrium fo(x,v,t), setting 
f(x,v,t) = f (x,v,t) + h 2 fi(x,v,t) in the semiclassical quantum Vlasov equation fl56l) and 



collecting equal powers of h . In other words set 

^>-^(5^)M^^)' (98) 

on the assumption of a non-degenerate quasi-equilibrium state. After calculating f\, insert 
f = fo + h 2 fi in the pressure in Eq. (165|) . A quantum hydrodynamical model similar to 
Eqs. flBBT) and (JHSHHZD is then found (and generalized by the inclusion of an energy trans- 
port equation). By definition, the resulting system is restricted to quasi-Maxwellian, dilute 
systems. We also observe that in the case of a self-consistent problem the electrostatic 
potential should have been also expanded in powers of some non-dimensional quantum pa- 



rameter, which was not done in 43j. The conclusion is that the procedures involving a 
Madelung decomposition of the quantum ensemble wavefunctions or a quantum corrected 
Wigner function equilibrium involve working hypotheses which are not rigorously justified. 
Starting from the Wigner-Poisson system, which is by definition collisionless, hardly one 
could derive rigorous macroscopic theories relying on quasi-equilibrium assumptions. Nev- 
ertheless, the numerical and analytical advantages of quantum fluid models over quantum 
kinetic models make them much more popular than the kinetic treatment, to approach the 
nonlinear aspects of quantum plasma physics. 

Alternatively, general thermodynamic arguments by Ancona and Tiersten have been used 
to derive the Bohm potential. Ref. 44| argues that the internal energy of the electron fluid 



in an electron-hole semiconductor should depend not only on the density but on the den- 
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sity gradient too, in order to extend the standard drift- diffusion model so as to include the 
quantum-mechanical behavior exhibited in strong inversion layers. Their method defines a 
"double-force" and a "double-pressure vector" which allows for changes of the internal ener- 
gy of the electron gas purely due to density fluctuations. Postulating a linear dependence of 
the double-pressure vector on density gradients (see Eq. (3.3) of [44j) and working out the 
conservation laws of charge, mass, linear momentum and energy, Ancona and Tiersten found 
a generalized chemical potential composed of two contributions: (a) a gradient-independent 
term which can be modeled by the equation of state of a zero-temperature Fermi gas or any 
other appropriated form. This corresponds to the classical pressure P c (n) of the quantum 
hydrodynamical model for plasmas; (b) a Bohm potential term proportional to a phenomeno- 
logical parameter left initially free. 



Later, Ancona and Iafrate 



of Ref. 



44j. In 



451 ] obtained the expression of the phenomenological coefficient 



451 ] , which is similar to 



43| , the first order quantum correction for a Maxwell- 
Boltzmann equilibrium was found from the semiclassical quantum Vlasov equation and then 
employed to calculate the particle density and the stress tensor. Eliminating the potential 
function between the two expressions, an equation of state relating the stress tensor and 
the particle density and gradient was derived, containing the Bohm potential. The gradient 
dependence of the stress tensor arises because of the quantum mechanical nonlocality. Once 



again, the demonstration in [45| is valid when quantum contributions to the self-consistent 
mean field potential can be ignored: no expansion of <fi in powers of a quantum parameter 
was used. Moreover, the dilute and semiclassical situation is supposed, where the Maxwell- 
Boltzmann statistics apply and a small dimensionless quantum parameter exist. 

Further, quantum hydrodynamic models with a "smoothed" potential have been derived, 



to handle discontinuities in potential barriers in semiconductors 46] . Energy transport 



and quantum drift-diffusion 



47] 



481 ] have been also employed to model transport in ultra-small 
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electronic devices, including the addition of viscosity effects 49]. 

It is amazing to notice how quantum hydrodynamic models are widespread in several 
areas, without apparently the explicit knowledge of the different communities. For instance, 
in molecular physics, the simplest gradient functional theory is given by the Thomas-Fermi- 
Dirac-von Weizsacker functional, which provides a suitable expression for the energy density 
of molecules 



51 



50j. It is composed by a Thomas- Fermi term for the uniform electron gas 



52| . a Dirac exchange functional 



53| . the Cou 



omb term for the electron g nuclear 



54|, with an adjustable fitting parameter. 



attraction term and the von Weizsacker term 
With this phenomenological parameter set to unity, the von Weissacker term gives exactly 
the Bohm potential. 

It is useful to remember the obvious limitations of the quantum hydrodynamic model 
for plasmas: (a) since it is a fluid model, it is applicable only for long wavelengths. This 
requirement is expressed roughly as A > A^> for non-degenerate and A > \f for degenerate 
systems, where A is the pertinent wavelength. Kinetic phenomena (Landau damping, the 
plasma echo {55)]) requiring a detailed knowledge of the equilibrium Wigner function certainly 
deserve a kinetic treatment; (b) no energy transport equation was included. This limitation 
is not definitive since it is a simple exercise to calculate it taking the second-order moment 
of the quantum Vlasov equation; (c) it applies only to non-relativistic phenomena; (d) no 
spin effects are included, except for the equation of state which can, in a certain measure, 
represent some quantum statistical effects. It suffices to apply the equation of state for 
a dense, degenerate electron gas. In the same way, some relativistic effects can be also 
incorporated with an equation of state for a relativistic electron gas; (e) no magnetic fields 
are in the present formulation; (f ) the sum of kinetic and osmotic pressures is assumed to be 
a function of the density only. This is in the spirit of density functional theories and appears 
reasonable in the vicinity of homogeneous equilibria, but may prove inappropriate in far 

31 



from equilibrium situations; (g) strongly coupled plasmas are out of reach since the starting 
point for the derivation was the quantum Vlasov equation, which is valid if correlations are 
not so important. 



37j starting from the Wigner- 



It is remarkable that the inclusion of magnetic fields 
Maxwell system can also be made in terms of the eikonal decomposition method, with the 
resulting quantum hydrodynamic model having just the addition of a magnetic ~ u x B 
force. 



V. CONCLUDING REMARKS 



Quantum plasma theory started in the 50's with the analysis of the ground state and 



correlation energies of the dense electron 
theoretical techniques. Pioneering works 



gas in metals 



17 



20 



56H58], through quantum field- 



59| also focused on collective oscillation 



modes in dense plasmas, described by the collisionless quantum Boltzmann equation. After 
a more or less quiescent period, in the last decade we observe a new wave of interest on 
quantum plasmas. Besides applications to diverse systems ranging from nanoscale electronic 
devices and dense astrophysics environments to intense laser-solid density plasma interaction 
experiments, an important step justifying the renewed attention to quantum plasmas is the 
emergence of efficient macroscopic models. This short introduction has been essentially 
dedicated to the analysis of the transition from the microscopic (kinetic) to macroscopic 
(fluid) modeling of quantum plasmas. 

Quantum plasma physics is a rapidly growing field, with unpredictable consequences. 
Below we list other fundamental topics for brevity not discussed here, currently being scru- 
tinized in the area. 

• The dynamics of spin degrees of freedom and the underlying ferromagnetic behavior 
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in quantum plasmas [60|-[64]. 



Relativistic quantum plasma models 



65]- 



681 ] . an essential subject since frequently 



quantum and relativistic effects appear simultaneously as e.g. in applications to com- 
pact astrophysical objects and intense laser plasmas. 

Extensions to include exchange-correlation and collisional contributions. If the anti- 
symmetric nature of the N— bod y w avefunction is considered a Hartree-Fock term 



appear in the interaction energy j20|. However, frequently such nonlocal exchange- 
correlation terms are replaced by local phenomeno logical expressions, along with the 



adiabatic local density approximation 



69 



70[. Also, we mention the Wigner-Fokker- 



Planck system approach to model quantum dissipation 



34 



Least but not last, in this introductory text we have focused on the detailed derivation 
of the basic kinetic and hydrodynamic models for quantum plasmas, in such a way to 
be accessible to beginners. A more complete account on the present applications of such 



quantum plasma models can be found in 



1] 



4]- 
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